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Abstract 

Although the neutron (n) does not carry a total electric charge, its charge and magnetization 
distributions represented in momentum space by the electromagnetic form factors, F^\q^) and 
F2^\q^), lead to an electromagnetic potential of the neutron. Using this fact, we calculate the 
electromagnetic corrections to the binding energy, i?^, of the deuteron and a one neutron halo 
nucleus (^^Be), by evaluating the neutron-proton and the neutron-charged core (^°Be) potential, 
respectively. The correction to Bd (~ 9 keV) is comparable to that arising due to the inclusion 
of the A-isobar component in the deuteron wave function. In the case of the more loosely bound 
halo nucleus, ^^Be, the correction is close to about 2 keV. 

PACS numbers: PACS: 21.10.Dr, 13.40.Ks, 13.40.Gp 
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I. INTRODUCTION 



Studying the sub-structures of nuclei and nucleons, encoded either in elastic electromag- 
netic form factors or deep inelastic structure functions, has been a field of continued interest 
in nuclear and particle physics. In this paper we shall focus on the aspect of nucleon sub- 
structure in connection with the appearance and relevance of electromagnetic form factors 
in nuclear physics. We shall calculate the electromagnetic contributions to the binding en- 
ergies of loosely bound neutron-nucleus systems, arising due to the neutron form factors. 
Such systems are realized by the deuteron and the one-neutron halo nuclei. Recently, the 
availability of radioactive beams opened up the possibility to study the structure of unstable 
nuclei. Such experiments Q revealed neutron rich nuclei whose spatial extensions are very 
large as compared to the range of the nuclear force. These so-called "halo" nuclei consist of 
very loosely bound valence neutrons which tunnel to distances far from the remaining set 
of nucleons. Such nuclei can be viewed as a system of a "core" with normal nuclear density 
and a low density halo of one or more neutrons. Thus one expects that the neutrons in the 
halo do not experience the strong force due to individual nucleons in the core but rather 
interact with the core as a whole. Based on this understanding, manyfew body models have 
been constructed and elaborately refined over the past few years 0,1^ in order to explain a 
variety of experimental data. 

It is well known by now that the neutron does have structure and hence an electric charge 
distribution which can be measured in elastic electron-nucleon scattering experiments. The 
interest in this field was revived due to the new experimental results on the nucleon form 

n 

factor from the Jefferson laboratory Considering the recent interest in nucleon form 
factors along with the experimental advances being made in halo nuclear studies, we found 
it timely to investigate the role of the neutron structure in the deuteron and one-neutron 
halo nuclei. In what follows, we shall derive an expression for the electromagnetic potential 
between the neutron and a charged particle (generalizing hereby the Breit equation by adding 
to it the finite size corrections) and then apply it to determine the perturbative corrections 
to the binding energy of nuclei. 
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II. THE NEUTRON ELECTROMAGNETIC POTENTIAL 



There exists a well known prescription for obtaining a potential from quantum field theory 
0,0]. The method consists of taking the Fourier transform of a non-relativistic scattering 
amplitude, say Mjyji, for a scattering process of the type AB AB. Since the method is 
completely general, one can handle simple cases with an amplitude which corresponds to one 
or two particle exchange diagrams or consider more complicated cases where the amplitude 
needs to be calculated using higher order corrections in the perturbative Feynman-Dyson 
expansion. If the amplitude depends only on the magnitude of momentum transfer {Q) in 
the process, then the Fourier transform to obtain the potential, V"(r), reduces to 7|: 

1 1 /""^ 

V{r) = ^-471 dQQ M^niQ^) sin(Qr) . (1) 

Obviously, if Mnr^Q"^) oc l/Q^ (non-relativistic propagator of a massless particle), the 
potential is proportional to l/r. Some examples of exotic potentials derived from quantum 
field theory can be found in ja fl Q • 

The aim of the present work is to study the role of the neutron structure in neutron - 
charged particle interactions. Hence we shall be interested in obtaining an electromagnetic 
potential which describes the interaction between the charged particle and the neutron. To 
obtain this potential, we shall start with the scattering amplitude for the process n + A ^ 
n + A, where A can be any charged particle with charge Ze (with Z a positive integer). We 
shall consider two cases: the neutron (spin 1/2) + A (spin 1/2) and neutron (spin 1/2) + A 
(spin 0) case with the former being relevant for the neutron-proton system (deuteron) and 
the latter for the one- neutron halo ^^Be taken as a neutron plus a ^°Be core in the s-state. 

Mathematically, the complete spin-1/2 fermion vertex (fermion-photon-fermion) is given 
as ^: 

2m / 4m f 

+ T^U' - 7^r)l5F,{q') , (2) 



2m/ V 2m/ 

where the sub-structure of the fermion is contained in the various form factors, Fi {i = 
1,2,3,4). is the four-momentum carried by the photon (g^ = q^q^)., rrif is the mass of 
the fermion and the 7's are the usual Dirac matrices jl^. The four-momentum squared. 
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= cu^ — (where uj is the energy and Q the three momentum), reduces in the non- 
relativistic hmit {uj ~ 0) to = —Q^. The above four form factors appear in the expression 
for the cross section for electron-nucleon elastic scattering and the nucleon form factors are 
thus extracted from such scattering experiments. The first form factor, Fi is connected 
to the charge distribution inside the nucleon and Fi{0) = the total charge of the nucleon. 
The form factors F2, F3 and F4 are related to the anomalous magnetic moment, electric 
dipole and Zeldovich anapole moment, respectively. In what follows, we shall obtain the 
electromagnetic potential for the spin-1/2 - spin-1/2 and the spin-1/2 - spin-0 case in terms 
of the nucleon electromagnetic form factors Fi{q'^) and F2{q^) as the other form factors 
{F^{q^) and F^{q^)) are related to the weak interaction and their effects are small. 

A. The neutron proton case 




FIG. 1: Feynman diagram for neutron scattering from a nucleus with Z being the number of 
protons. 

In order to perform a complete calculation of the potential, we turn to the task of expand- 
ing the amplitude in 1/c^ terms, thereby generalizing the Breit equation 5] by the inclusion 
of two electromagnetic form factors. This will lead to non-local and spin-dependent terms in 
the potential whose contributions to the binding energies are not negligible as will be seen 
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later. Starting with the Feynman diagram of Fig. 1, the nucleon-photon-nucleon vertices in 
terms of the nucleon form factors Fi and F2 can be written as, 



cr 



2Mxc 



X 

2 



(3) 



= F^Y + 



a 



QuF^. 



2Myc 

The photon momentum, q = P[ — Pi = P2 — P-^. In the non-relativistic hmit {q^ = 0), 
= —Q^, where, Q = — Pi = P2 — P'z- The amphtude for the process N + N —>■ N + N 
is then given by [5|, 



M 



fi = ex ey 



(4) 



where, D^y{Q) is the photon propagator and u{pi), u{p2) etc., the usual Dirac spinors given 
as, 



u 



1 



p' 



w 



a ■ p 

2Mc 



w 



(5) 



Substituting for the u{p)^s and the vertex factors given above, the amplitude is evaluated 
and then rearranged to be written in the following form: 



Mfi = -2Mx ■2My [w'* w'2* Uipi,p2, Q) wi W2 \ 



(6) 



thus obtaining the potential U{pi,p2,Q) in momentum space. The potential in r-space is 
obtained simply from a Fourier transform of the potential in momentum space, namely, 

d^q 



V{r} 



e^^-^ U{pi,p2,Q) 



(7) 



(2 7r)3 ■ 

The electromagnetic neutron-proton potential in momentum space, using the above proce- 
dure is found to be, 

1 1 1 (Ti ■ (72 

^ ~ 8M2c2 ~ 8 M2 c2 " 4M„Mpc2 

iai ■ {Q X P2) 



^ • Q)i^2 ■ Q) _^ iffl • {Q X Pi) _ 02 ■ {Q X P2) _^ 02 ■ {Q X pi] 



4 Mn Mp d^Q^ 



+ F^Fi 



jpi ■ Q) iP2 ■ Q) 

(?1 • 0^2 



^ i^l ■ Q){^2 ■ Q) ^ ■ {Q X Pi) 02 ■ {Q X P2) 



AM^c^ AM^MpC^ AMnMpC^Q^ 2MpMnC^Q^ 2M^c^Q^ 
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+ Fi 



01-82 ^ {8i-Q){a2-Q) 



AMnMpC^ AMnMpC^Q^ 

In obtaining the above expression, we have dropped terms of order higher than (1/c^). 
The above expression involves a separate 1/c^ expansion corresponding to each of the form 
factor combinations. This is to say, the term proportional to F^Ff /Q'^ is the leading term 
of the coupling constants combination F^Ff, whereas the other terms proportional to this 
combination are suppressed. However, we cannot say that just like the (1/c^) terms in 
F^Ff, those in the other combinations like -F^-Ff') ^^"^ F2P2 ^-^'s also suppressed, as 

the latter represent different combinations of coupling constants. Indeed, all combinations 
different from F"Ff are displayed also in the leading order albeit they appear with the 
factor which does not imply that they are smaller than F^FljQ'^. We will see later how 
this works in practice. The constant e in the above is defined such that = a, the fine 
structure constant. Since the experimental form factors for the neutron as well as the proton 



are normalized in a similar way, we have a vertex factor of v47ra at the neutron-photon- 
neutron vertex too, in order to remove the extra factor by which the experimental neutron 
form factor was divided. 

The detailed calculation of the Fourier transforms of the terms in (jH)) and the corre- 
sponding potential in r-space is given in the appendix. The terms which appear with cross 
products of momenta in (jH}, do not contribute to the binding energy correction for the 
deuteron, as shown in the appendix. 



B. The neutron spin-zero-nucleus case 

To evaluate the electromagnetic potential between the neutron and a spin-zero nucleus, 
we repeat a very similar procedure as above, with the vertex Y in Fig. 1 replaced in this 

case by, 

r^ = F^(g2)(P2 + P^r, (9) 

where, F^(g^) is the form factor of the nucleus (A) in momentum space. We also replace 
the appropriate normalization for the spin zero particles. The potential in momentum space 
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is found to be, 



f/"^(pi,P2,Q)=47re^^ F^F 



1 



Pi ■P2 



+ 



iai ■ (Q X pi) 



+ 



(Pi • Q) {P2 ■ Q) 
iai ■ {Q X P2) 



+ F^F 



+ 



4M2c2g2 2MAMnC^Q'^ 
idi ■ {Q X pi) ■ {Q X 



(10) 



4M2c2 2M2c2g2 2M^M„c2Q2_ 

This potential will be used in the present work to evaluate the correction to the binding 
energy of the one-neutron halo nucleus, ^^Be, taken as a neutron plus ^°Be core. The 
evaluation of the ^^Be wave function is usually done in a neutron plus ^°Be core model 
where the ground state (spin zero) and an excited state (spin 2) of the core are considered. 
We shall restrict to assuming the core to be a spin zero nucleus since the calculation of a 
spin 1/2 - spin 2 electromagnetic potential is beyond the scope of the present work. The 
contributions of the spin-2 terms to the ^^Be wave function are in any case small and would 
have only a small overlap with the electromagnetic potential which is not long-ranged. 

In the next sub-section, we shall describe the various parameterizations of the form factors 
which will be used to evaluate the above potentials. 



C. The electromagnetic form factors Fi{q^) and F2{q^) 

The form factors Fi and F2 have been extracted from several experiments and parame- 
terized in different forms in the past However, a discrepancy between the old results 
and new experiments which extract the two form factors through polarization measurements 
was recently reported . It was found, however, that this discrepancy could be resolved by 
taking into account the two photon contribution [l^. Therefore, in the following, we shall 
make use of the 'standard' parameterization of the form factors. Considering the renewed 
debate on the nucleon form factors, in the present work we shall also perform the calcula- 
tions with different parameterizations available in the literature. The Sachs form factors Ge 
and Gm which appear in the expressions for the electron- nucleon elastic cross sections are 
related to the structure functions Fi and F2. The nucleon form factors F^ {q^) and F^{q^) 
are thus given as 

r,N. 2^ AM'G^{q')-q'GUq') .... 
^1 ) = 4M2 - g2 ' (11) 
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4M2 - g2 

where M is the mass of the nucleon and = uj"^ — is the four momentum of the 
virtual photon as mentioned before. A large body of experiments starting from the sixties 
until now has been dedicated to the extraction of G'^(g^) and G\j{q'^). Though the newer 
experiments have smaller statistical errors, there still exist uncertainties arising from the 
theoretical description of the deuteron. Hence, in the case of the neutron, we shall use 
different parameterizations available in literature. 



1. The dipole form factor and the neutron electric potential 



We shall use the following form for G'%{q^) [13j with different sets of parameters 'a' and 
'6' obtained in literature: 



(12) 



1 - 6(gV4M2) 

where Hn is the neutron magnetic moment and Goiq"^) is the standard "dipole fit" which 
is generally used in summarizing the electron-nucleon elastic scattering data. With G£,{q^) 
defined as 



[1 — q^/m?)"^ 



it was observed that 



GWl^G^^ 



GD{q' 



(13) 



(14) 



where the magnetic moments fip and of the proton and neutron (in nuclear magneton) 
are 2.79 and —1.91, respectively. Using the above definitions of (j^(g^) and (j^^(g^) in (fTTjl . 
we obtain the following non-relativistic expression (i.e. with q^ ~ —Q"^) for the form factor 
Ff(Q2): 

"4M2 (a - 1) - Q%' 



1 



1 



4M^ + Q^ (1 + QVm- 



2^2 



(15) 



4M2 + 6Q2 

At this point, it is nice to note that using the dipole form factors, the potential between 
the structured neutron and an electron can be derived analytically and is given by, 

-Za 1 



y"(r) 



(1-6) r 



{a-l)H^ -hH'^ 



(16) 
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where 



.12 



+ 



.12 



and 



1 



- 1) 
1 



^—Kmr _j_ 



f2 



r e 



.12 



K 



'\ 2 



12 



1)^ 



1)^ 



1 



1 



f«:2 - r 



,/2 



1 



(17) 



(18) 



/€2 V 2 

where k = 2M/m, k' = n/Vb and typically = 0.71 GeV^ HQ. 
In the neutron magnetic form factor was also fitted with the form 

Gi,m 1 

Hn 1 - l-74g + 9.29g2 _ 7.63g3 + 4.63Q4 ' 
and was shown to reproduce the existing data quite well. We shall perform calculations using 
this form of G'%j{Q'^) and the dipole form of G'%{Q'^). In this case however, we evaluate the 
potential numerically and use it to evaluate the binding energy as explained in the next 
section. 



(19) 



2. The hump tail model 



All new and old data of the electric and magnetic form factors of the proton and the 

n 

neutron were recently re-analyzed by a fit [1^ based on a pion cloud model. A common 
feature (a bump) was noticed in the data at low momentum transfer, which the authors 
attributed to the existence of a pion cloud around a bare nucleon. The form factors were 
parameterized in terms of a bump on top of a large smooth part. Purely phenomeno logically, 
the smooth part was parameterized with a superposition of two dipoles: 



+ 



(l + QVan)2 (l + QV«2i)2 
and the bump as a superposition of two Gaussians: 



exp 



1 ( Q-Qb 

2 \ ab 



+ exp 



1 ( Q + Qb 

'2 I ab 



(20) 



(21) 



The neutron and proton, electric and magnetic form factors are then given by the ansatz, 

Gn{Q^) = GsiQ') + a, Q2 g,{Q'') (22) 



where the parameter is essentially the amplitude of the bump. The parameters for the 
proton form factors are taken from Table II of ^ and those for the neutron are taken from 
a more recent work 16) where the authors use a similar phenomenological fit. The neutron 
parameters fitted in (l^ are within the error bars of those quoted in Table II of for the 
neutron. 



D. The nuclear form factor 

The form factor of ^°Be is evaluated by taking the Fourier transform of a semiphenomeno- 
logical nuclear charge density given in (l7|. Though the actual density of ^°Be may have a 
complicated structure, the following form provided by the authors in is in very good 
agreement with experimental data on nuclear form factors. The charge density distribution 
is given as, 

^(^) = 1 + [(1 + (r/i?)2)/2]"[e(^-^)/'^ + e-(^'+^)/«] ' ^^^^ 
where, p° is determined from the normalization condition: 

47r ( p{r) r^dr = Z , (24) 



with Z being the total 



found in Table 1 of 



17|. 



number of protons in the nucleus. The parameters for ^°Be can be 



III. CORRECTIONS TO NUCLEAR BINDING ENERGIES 

As mentioned in the previous section, we evaluate the correction AE to the binding 
energy of the deuteron and the one-neutron halo nucleus ^^Be. 

AE = y ^*(r)V"(r)^(r)dr, (25) 

where \1/ is the wave function corresponding to the unperturbed Hamiltonian, Hq (the total 
H = Ho + V). 
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A. The deuteron 



In the deuteron case, we perform the calculation for the dominant s-state of the wave 
function obtained using the neutron-proton strong interaction. When the potential, V"'{r), 
depends only on the magnitude of r and is spin-independent, the above equation reduces to 
the simple form: 

roc 

AE = r^u^{r)V'''{r)dr (26) 



where u{r) is the radial part of the deuteron wave function. We use a parameterization 
of this wave function using the Paris nucleon-nucleon potential as given in il8j. Since the 
Paris potential itself is written as a discrete superposition of Yukawa type terms, the wave 
function is parameterized in a similar way as: 

13 

u{r) = ^ Cj exp(— mj-r) / r . (27) 



The coefficients Cj with dimensions of [fm~^/^] are listed in Table 1 of 18 1 and the masses 
rrij are given as rnj = a + {j — l)mo, with niQ = 1 fm~^ and a = 0.2316 fm^^. In Table I, 
we list the corrections to the deuteron binding energy due to each of the terms in (jH)) using 



the most recent parameterization of Refs. 



and |l6|. It can be seen that with the F^Ff 



combination of form factors, all terms other than the leading one are one or two orders of 
magnitude smaller as expected from the (1/c^) expansion. The contributions of the terms 
involving F^((5^) are large as compared to those with F"(Q^) simply due to the fact that 
unlike F"((5^), -F^((5^) does not vanish as 0. Table I has been arranged according to 

the different form factor combinations to display in a clear way, the theoretical expectations 
discussed below Eq. (jS)). First, one should note that the entry corresponding to the F"Ff 
form factor combination has a leading term followed by (1/c^) suppressed next-to-leading 
order terms. We have done this to display in this particular case, the significance of the 
next-to-leading order terms. In the next three combinations, only the leading terms of these 
particular combinations are listed. These terms contain a (1/c^) factor which, however, does 
not make them always suppressed as compared to the leading term of the first combination 
FJ^Ff, bearing in mind that we introduce a new momentum dependent coupling constant 
F2. This is to say that one should not expect terms with different form factor combinations 
to give similar results. This is in agreement with the remarks we made regarding the 1/c 
expansion in section III Al 
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TABLE I: Corrections due to individual terms in the neutron-proton electromagnetic potential, to 
the deuteron binding energy in keV, using the parameterizations of form factors in [isl . Il^ . The 
corrections assuming a point proton are listed in the last column. 



Form factor 


Term in U^^{pi^p2-,Q) 


AF (keV) 


AF (keV)(with Ff = 1, F| = Kp) 


combination 








Ff Ff 


1 


-1.25 


-0.923 




1 1 


-0.016 


-0.049 




8M2c2 8M2c2 




0\ ■ (72 

4 M„ Mp c2 


-0.016 


-0.049 




(<Ji • q){a2 ■ q) 


-0.028 


-0.063 




Pi 

MnMpC^q^ 


-0.02 


0.053 




{pi ■Q){P2-q) 

MnMpC^q^ 


-0.05 


-0.155 


F^F^ 


1 

4M2c2 


-0.034 


-0.177 




(71 • (72 

4 M„ Mp c2 


-0.034 


-0.177 




(ai • q){a2 ■ q) 
AMnMpc'^q^ 


-0.039 


-0.1135 


F^Ff 


1 

4M2c2 


3.54 


2.986 




ai ■ (T2 

4 Mn Mp c2 


3.54 


2.986 




(ai • q){a2 ■ q) 
AMnMpC^q^ 


-1.32 


-2.273 


F^FP 


d\ ■ (72 

4 Mn Mp c2 


6.56 


5.355 




(<Ti • ^((72 • q) 
AMnMpC^q^ 


-2.02 


-4.076 




Total: 


8.81 


3.3735 



In the last column of Table I, we list the corrections assuming the structured neutron 
and point proton interaction. The potential in this case increases a lot in magnitude (for 
example, the depth of the potential in the leading term is about —270 keV assuming a point 
like proton as compared to the —60 keV in the structured proton case), however, the overlap 
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of the deuteron wave function and the potential does not change drastically. 

In order to demonstrate the typical form of the neutron-proton electromagnetic potential, 
in Fig. 2, we plot V"'{r) as a function of r for two terms which contribute to the binding en- 
ergy correction with opposite signs. The solid line corresponds to the attractive potential of 
the Ff (Q^) leading term and the dashed line to the spin-independent F2{Q'^) Ff^Q"^) 

term. The other major contributors to the total correction /S.E are spin-dependent and con- 
tain operators. However, if one evaluates the effective potentials corresponding to these 
terms as mentioned in the appendix, these potentials are seen to have a similar form and 
range as the ones shown in the figure. The effective potential corresponding to the biggest 
contribution coming from the F^(Q^) ^^'(Q^) is repulsive with \/"(r — >^ 0) ~ 220 keV. 





Terms containing 



f;(q') f;(qVq' 
fAq') f;(qV(4m;c') 



[b] 



2 

r (fm) 



FIG. 2: (a) The I = radial wave function squared of the deuteron using the Paris potential 
and (b) two different terms of the electrom agne tic potential between the neutron and the proton 



evaluated using the parameterization in 
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In Table II, we list the total Ai? due to all terms using the parameterization of Ref. 



with different choices of the parameters a and h in 
with a = corresponds to the option G^{Q'^) = 



). The first parameter set in this table, 
1911. The second set with a = 1, 6 = 5.6 
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TABLE II: Corrections to the deuteron binding energy, AE in keV, using the parameterization of 
form factors as in Q. m(i'^) taken as in Eqs (4), (5) and G%[{q^) as in Eq. (7) in Q. The 
parameters a and h appear in G^{Q'^). 



a 


b 


AE (keV) 
(Total) 







7.9 


1 


5.6 


11.24 


1.25 


18.3 


11.8 


1.12 


21.7 


11.39 


0.94 


10.4 


11.01 


0.94 


11 


11 


0.9 


1.75 


10.83 



is the well-known Galster parameterization |20|. The values in the next two rows lie within 
error bars of the best fit values obtained in [l3|, namely, a = 1.25 ± 0.13 and b = 18.3 ± 3.4. 
The next two rows lie within the error bars of a = 0.94 and b = 10.4±0.6 which were 
obtained by constraining the slope of to match the thermal neutron data Q]. Finally, 
the last value corresponds to the one used in in connection with the electromagnetic 
N —>■ A transition. 

It is interesting to note that these corrections to the deuteron binding energy are com- 
parable in order of magnitude to those due to the presence of a A-isobar component in the 
deuteron. In 2^ , the corrections due to the A were found to be around 3 keV. The deuteron 
is a precision tool of nuclear physics, both experimentally (the binding energy is precisely 
known and theoretically (as it is a two body problem). The level of accuracyidemands 
even to discuss and agree upon corrections of the order of eV. For example, in [2^], the so 
called "Doppler broadening of the 7 ray" arising basically due to the kinetic energy of the 
neutron in the ^H(n,7)^H reaction (used to extract the deuteron binding energy) is found to 
introduce an error of 25 eV, maximally. This correction is considered by the authors in j^l , 
to be a large one as compared to the 2.3 eV error which comes from the 7-ray detector itself. 
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The 10 keV electromagnetic correction is then quite large as compared to the above correc- 
tions. To obtain more accurate estimates of the electromagnetic corrections, it is therefore 
clear that there is a need to know the neutron form factors to a better precision and for a 
larger range of the momentum transfers than what we know today. 

There is a non-zero, albeit small effect of this electromagnetic correction to the physics of 
nucleosynthesis. Assuming that some fundamental constants, among them the fine structure 
constant a, change with the cosmological epoch 2^, this variation can affect the abundances 



of primordial light nuclei. In [26[ such an analysis for the deuteron gives 

— ^ = 2.320— + ... (28) 

Yd a 

where the dots indicate the contributions of the variations of other constants. The numerical 
value 2.320 above includes all sources of the a dependence of the deuteron abundance, 
except for the direct dependence of the binding energy B^, on a. To estimate this effect on 
6Yd/Yii we make use of the equilibrium solution for the deuteron abundance given by 
Yd oc expl—Bd/ksT] and write Bd = + CiB'^m = -^d + -Bem, to obtain 

' 5Yd\ 5a ( B^rn\ j^29) 



"^Y d J fWvc^r'^ ^ ''^ 



direct 

As an estimate we take, ksT ~ 5 x 10^^ MeV, corresponding to the point where the 
equilibrium solution deviates from the numerical one due to the depletion of the deuteron 
in the production of heavier nuclei. Note however, that the final abundance does not differ 
much from the equilibrium solution. This gives Bd/ksT ~ 0.2 which changes the coefficient 
2.320 in to 2.120. The analysis of 6Yd/Yd helps in searching for positive signals of the 
variation of fundamental constants. 



B. One neutron halo - Be 

The wave function for ^^Be is taken from a coupled channel calculation performed for 
one neutron halo nuclei. A deformed Woods-Saxon potential for the neutron-core interaction 
is used to take into account the excitation of the ^°Be core. The coupling of the neutron 
to the i°Be core (taken as 0'^ or 2+) gives rise to the three components, 2si/2, 1(^3/2 and 
ld5/2 of the ^^Be wave function. The notation used is nlj with / and j being the orbital and 
total angular momenta, respectively. The normalization of the radial wave functions is such 
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TABLE III: Corrections to the ^^Be binding energy, Ai? in keV, using the parameterization of 

[l4|. Gl^{q^) is as in Eq. (7) in Q- 



neutron form factors as in 



G%{q ) is as in Eq. (7) in [ij]. The parameters a and h appear in 



a 


h 


(keV) 


AS (keV) 


AS (keV) (Total) 






due to the 


due to the 








term containing 


term containing 








Ff(Q2) 











-3.52 


3.53 


0.01 


1 


5.6 


-0.54 


3.55 


3.01 


1.25 


18.3 


-0.12 


3.56 


3.44 


1.12 


21.7 


-0.54 


3.55 


3.01 


0.94 


10.4 


-0.82 


3.55 


2.73 


0.94 


11 


-0.83 


3.55 


2.72 


0.9 


1.75 


-0.74 


3.55 


2.81 




Ref.[l^,i^ 


-1.32 


3.17 


1.85 




Ref.[15,i6] (FioBe = 4) 


-3.93 


17.93 
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that j[u{r)fdr = 0.85,0.02 and 0.13 for the 2si/2, 14/2 and 1^5/2 waves, respectively. The 
total wave function for ^^Be as given in is expressed as a sum over spins written in terms 
of the rotational matrices. However, since we have evaluated the potential for the spin-1/2 
spin case of a neutron and nucleus, we shall be presenting results only due to the s-wave 
component. All components of the wave function are however plotted in Fig. 3, where one 
can see that at least for the leading spin-independent term, the d-waves would have a much 
smaller overlap with the potential as compared to the s-wave. 

For a detailed comparison of this neutron-^°Be core model with the shell model for ^^Be 
as well as experiments, we refer the reader to the work of F. M. Nunes et a/ j^. 

In Table III, we list the corrections to the ^^Be binding energy of 500 keV using different 
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parameterizations of the nucleon form factors as in the case of the deuteron. We also hst the 



results in the parameterization of 



15 



, 1^ for the case assuming a point like charged nucleus. 
One can see that the binding energy correction turns out to be an order of magnitude larger 
than the one which takes the nuclear structure into account. This big difference (as compared 
to the not so large one in the deuteron case with point like proton) can be understood in 
terms of the fact that the nuclear form factor falls more rapidly with momentum and the wave 
function of the halo nucleus has a different behavior as compared to that of the deuteron. 
Thus, if the nuclear form factor is not included, even the order of magnitude of the resulting 
corrections is incorrect. 

Note that in performing the calculation of the electromagnetic correction for the one- 
neutron halo, we assumed the same approach as taken for the strong interaction part, i.e., 
we treat the neutron to be interacting with the rest of the nucleus as a whole. 

0.2 




> 



> 



30 
20 

10 - 

-10 



Terms containing 

— f;(qV,„^^(qVq' 
^\-- f;(q')f,o (qV(4m;c') 



[b] 







2 4 
r(fm) 



FIG. 3: (a) The radial wave function components squared for ^^Be within the neutron-^'^Be-core 
model 12]. (b) the electromagnetic potential between the neutron and the ^''Be core, evaluated 



using the parameterization in 
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IV. SUMMARY 



To summarize the findings of tlie present work, we can say tliat we liave adopted a new 
approacli to tlie problem of calculating electromagnetic corrections to the binding energies 
of loosely bound neutron systems encountered in the deuteron and one-neutron halo nuclei. 
The relevance of such corrections for the deuteron lies in the realm of precision nuclear 
physics as -B^ is an accurately known number. A small correction to the analysis of the 
variation of primordial deuteron abundances with respect to the change of fundamental 
constants can also be found. To be able to perform the calculation with a reasonable 
accuracy we generalized the Breit equation by allowing the couplings constants to vary 
with momentum transfer (form factors). To compare the result with other loosely bound 
systems, we calculated the electromagnetic correction for the one-neutron halo nucleus too. 
This calculation derives its relevance also from the exotic nature of the halo nuclei. Given 
the low binding energies of the one-neutron halos and the fact that in the standard picture, 
the valence neutron resides far from the remaining core, one could naturally ask the question, 
does the electromagnetic interaction play some role apart from the strong interaction which 
is not so strong as in normal nuclei? The answer is that the electromagnetic corrections 
viewed as a ratio to the corresponding binding energies are of the same order of magnitude 
in the case of the deuteron and the one-neutron halo provided we treat the latter as a neutron 
interacting with the rest of the nucleus (we mentioned already that this is in the same spirit 
as the calculation done for the strong interaction part). The strength of the electromagnetic 
potential between the neutron and the core depends on the number of protons in the core. 
Hence, such a correction in the case of the one-neutron halo ^^C would be larger; more so 
due to the fact that the binding energy of the ^^C can be smaller ~ 240 MeV |2q . 

All this also shows that we ought to know the electromagnetic form factors of the neutron 
more precisely. 
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APPENDIX A: THE POTENTIAL IN COORDINATE SPACE 



As mentioned in the main text, the potential in r-space is evaluated as a Fourier transform 
(FT) of U{'pi,'p2,Q). In what follows, we shall enumerate the terms in the order in which 
they occur in (jSI). The FT of the first three spin- independent terms in (|S)) is trivial and 
gives, 

V,{r) = ^ r F-{Q') dQ (Al) 

nr Jo Q 

V^2(3)(r) = — J / F^iQ') Ff Q sm{Qr) dQ . 

vrr V 8 M^(p) / Jo 

— * 

The fourth term which is spin-dependent is evaluated using the relation, S = (1/2) (ai + 0^2), 
leading to <?! • = 2 — 3. Hence, 

2eV 1 \ /-^ 

The spin operator S"^ operates on the deuteron spin function when one evaluates the binding 
energy correction as in (j2S)). Since the deuteron spin is, S* = 1, (2 5^ — 3) Xd = +1 Xd and the 
potential V4 gives a contribution similar to V2,3 up to a factor of about 2. The evaluation of 
V^lr) involves a small mathematical trick as described in |5j]. The FT of this term is written 
as, 

V,{r) = J e^'^-^fiQ') ^ (ai • Q) {a, • Q) ^ , (A3) 
where f{Q^) = F^{Q^) Ff{Q^) e"^ / {AMnMpC^). Now, 



K.(r) ^ ^ (-^j;^) f Q sHQr) dQ (2 - 3) . (A2) 



^(r) = -zai- VF(r), (A4) 

with 



2, 47r d^Q 



Fir) = I e^^- fm (^2 ■ Q) ^ . (A5) 



Repeating the same trick again, we can write V^{r) as. 



V,{r) = -{a,. V) {a, ■ V) / e^<^-^^ % ^ . (A6) 



g2 (27r) 



If we use the relation, 2{S ■ V)^ — = (cti ■ V) ((T2 • V), and we choose, f = r z, then 
the potential involves the operator S'^, which in the case of the deuteron again acts on the 
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spin-1 wave function. We can then write an effective potential (which can be used as in Eq. 
for evaluating AE) as, 



J. J. 2 r'^ 

V,'"{r) = - dQ f{Q^ 

OTT JO 



sin((yr) ( -— + 



— cos(Qr) 



(AT) 



The next four terms which involve cross products do not contribute to the binding energy 
correction. For example, writing the sixth term as, 

and defining. 



(A8) 



4Q2 M2 c2 (2 



V6(r) can be expressed as. 



\/6(r) = -47re'[V/(rO-(aixpi) 



where, 



Writing 



-* d X ' 

V/(r) ■ ((Ti X pi) = g:^f{r) eijk Oij pik 



(A9) 



(AlO) 



(All) 



Pik oc 



d d dr Xk d 



dxk dr dxk r dr ' 
and the deuteron wave function with its spatial and spin parts as, \1/ = R{r) Yqq{ 
energy correction as in ()25p for s-waves becomes. 



X, the 



(A12) 



The factor pi ■ p2 in the next term is given in the center of mass system of the two nucleons 
p 



as, pi-p2 
Thus, 



'2 = h\'^, and acts on the deuteron wave function in the evaluation of AE. 



Vio{r) 



2e 



1 



sin(Qr) 



92 



M M 2 , ^liQ') FfiQ') """tt^ dQ 

IT r Mn Mp & Jo Q or'' 



(A13) 



The last of the Fi"(Q2) Ff (Q2) terms is evaluated using a similar mathematical trick as was 
used for V^lr). Here again, working in the center of mass system, 

Fi"(Q2)Ff(g2)e2] 47r d^Q 



Vnir) = -(pi . V) (p2 ■ V) / e 



Mn Mp C2 Q2 



(27r)= 



(p-V) 



2 / ^iQr 



Mn Mp C2 g2 



47r ci^g 



\3 • 



(A14) 
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Here, the momentum operator p = —iTiV acts on the wave function [5| while evaluating /S.E. 
The remaining terms in (jH} differ only in the form factor combinations and can be evaluated 
as above. 
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